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We calculate the complete next-to-leading order QCD corrections to the three-body 
decay of gluino into top-pair associated with a lightest neutralino in the minimal su- 
persymmetric standard model. We obtain that the LO and NLO QCD corrected decay 
widths of g -> tiy°i at the SPS 6 benchmark point are 0.1490 GeV and 0.1069 GeV 
respectively, and the relative correction is —28.2%. We investigate the dependence of 
the QCD correction to g — > tixi on tan j3 and the masses of gluino, scalar top quarks 
and the lightest neutralino around the SPS6 benchmark point, separately. We find 
that the NLO QCD corrections suppress the LO decay width, and the absolute relative 
correction can exceed 30% in some parameter space. Therefore, the QCD corrections 
to the three-body decay g — > tixi should be taken into account for the precise exper- 
imental measurement at future colliders. Moreover, we study the distributions of the 
top-pair invariant mass and the missing energy (E mtss ), and find that the line 

shapes of the LO distributions of M t i and E miss are not obviously distorted by the 
NLO QCD corrections. 
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I Introduction 



The supersymmetric theories [T] are the most motivated extensions of the standard model 
(SM). They predict that the SM particles have their corresponding superpartners (sparti- 
cles), and the CERN Large Hadron Collider (LHC) and the International Linear Collider 
(ILC) might provide the experimental facilities to confirm the existence of these new parti- 
cles. Among all the supersymmetric models, the minimal supersymmetric standard model 
(MSSM) with i?-parity conservation is the most interesting one which is well studied in recent 
years pQ. In the i?-conserving MSSM, there exists a stable lightest neutral supersymmetric 
particle (LSP) which is called the lightest neutralino denoted as Xv 

Supersymmetry (SUSY) must be broken in the practical world and the sparticle mass 
spectrum depends on the SUSY breaking mechanism. The fundamental MSSM parameters 
need to be determined from the precise measurement of the masses, production cross sections 
and decay widths of these superpartners. With these parameters we can reconstruct the 
SUSY breaking mechanism and probe the MSSM. 

Among all the sparticles in the MSSM, the two colored scalar quark (squark) chiral 
eigenstates qi and qn are the corresponding superpartners of the chiral quarks appearing 
in the SM. The physical mass eigenstates q% and q~2 are the mixtures of these two chiral 
eigenstates. The scalar partners of top quarks are expected to be the lightest squarks in 
supersymmetric theories. As the colored supersymmetric particles, squarks and gluino may 
be produced copiously in hadronic collisions. After the (pair) production of these particles 
they are expected to decay via cascade decay to lighter SUSY particles associated with quarks 
and/or leptons. Gluino may have many decay modes in the MSSM depending on the mass 
spectrum of SUSY particles [2j[3]. In general, gluino has three major decay channels, i.e., (1) 
via an off-shell squark to two quarks and a LSP, g — > qq + Xi, (2) via an off-shell squark to 
two quarks and a chargino or heavier neutralino, e.g., g — > q<? + Xi> an d (3) through a quark- 
squark loop to a gluon and a LSP, g — > g + xl- The lowest order decay width of the decay 
process g — > qqxl was previously evaluated in Refs.[2j HI [5]. In Ref.[l] the authors studied 
also the two-body decays g — > gx° n including the one-loop QCD corrections by adopting 



the leading logarithmic approximation. Recently, a complete study for all two-body decay 
channels of the gluino in the complex MSSM with full one-loop electroweak effects has been 
presented in Ref.jBjj. If t\ is relatively lighter than other squarks and the gluino is heavier 
(rrig > mf ), the decay channel g — > tixi could have major contribution to the total decay 
width of gluino [HE]. Therefore, the accurate calculations including the NLO corrections to 
the three-body decay g — >■ tixi are necessary. 

Currently, ATLAS and CMS experiments [8j [9] with 4.7 fb~ l data have severely con- 
strained the masses of the strongly interacting SUSY particles — squarks and gluinos. How- 
ever, in some SUSY scenarios these SUSY particles may still be rather light. One interesting 
possibility is the light stop scenario. In Ref.[T0] the authors propose that the recently devel- 
oped techniques for tagging top jets can be used to boost sensitivity of the LHC searches for 
the SUSY scenario in which the third generation squarks are significantly lighter than those 
of the first two generations. 

In this paper, we focus on the NLO QCD corrections to the three-body decay g — > fix® ■ 
In section II we present the tree-level calculations for the decay process g — > ttXi ■ In section 
III we provide descriptions of the analytical calculations of the NLO QCD corrections. In 
section IV we present the numerical results around the scenario SPS6 point as proposed in 
the SPA project [HI [12], and discuss the dependence of the QCD correction on tan/3 and the 
masses of gluino, top squarks and the lightest neutralino. Finally, we give a short summary. 

II LO calculations for g — > tixi 

In this paper, we denote the decay process g — > tixi as 

~g{px) ^t(p 2 )+i(p 3 )+xi(p4), (2-i) 

where pi, p 2 , an d Pa ar e the four- momenta of the gluino and the decay products, re- 
spectively. The leading order (LO) Feynman diagrams for this decay are displayed in FigJU 
where t s with the lower index s running from 1 to 2 represent the two top squarks t\ and t 2 - 
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Figure 1: The tree-level Feynman diagrams for the decay g — > ttXi ■ 

The LO decay width of the process g — > tixi can be obtained by using the following 
formula: 



The summation is taken over the spins and colors of initial and final states, and the bar 
over the summation recalls averaging over the spin and color of initial gluino. (i$3 is the 
three-body phase space element defined as 



After doing the integration in Eq. ( 12.21) over the phase space, we can get the expression of 
the tree-level decay width of the decay process g — > tixi ■ 

In the LO calculations, the intermediate top squarks are potentially resonant. We adopt 
the complex mass scheme (CMS) [13] to deal with the top squark resonance effect. In 
the CMS, the complex masses of the unstable top squarks should be taken everywhere 
in both the LO and NLO calculations. Then the gauge invariance is preserved and the 
singularities of propagators for real p 2 are avoided. The relevant complex masses are defined 
as /i~ = m~ — imiTi., where mf. (i = 1,2) are the conventional pole masses, F^. (i = 1,2) 
represent the corresponding total widths of the top squarks, and the poles of propagators 
are located at uf on the complex p 2 -plane. Since the unstable particles are involved in the 
loops for the 0(a s ) QCD corrections, we shall meet the calculations of iV-point integrals 
with complex masses. 




(2.2) 
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Ill NLO QCD corrections to g -> ttjft 

The calculations of the decay process g — > tixi in the MSSM are carried out in t'Hooft- 
Feynman gauge. In the QCD NLO calculations, we use the dimensional regularization (DR) 
scheme to isolate the ultraviolet (UV) and infrared (IR) singularities. The Feynman diagrams 
and the relevant amplitudes are generated by using FeynArts3.4 [14], and the Feynman 
amplitudes are subsequently reduced by FormCalc5.4 [15] . The phase space integration is 
implemented by using the Monte Carlo technique. The NLO QCD corrections to the decay 
process g — > ttx® , denoted as ATnlo, can be divided into two parts: the virtual correction 
from one- loop diagrams (AT virtua i) and the real gluon emission correction (Ar rea /), i.e., 

Tatlo = Tlo + AT nlo = Flo + AT virtua i + Ar rea ;. (3.1) 

Both the virtual and the real corrections contain IR singularities. These IR singularities 
exactly vanish after combining the virtual correction with the real gluon emission correction 
together. Then the decay width of g — > tix\ including the NLO QCD corrections is IR-fmite. 

III.l Virtual corrections 

A. Definitions of counterterms 

The renormalized virtual NLO QCD corrections to the decay g — > tixi in the MSSM 
include the contributions from self-energy, vertex, box, and the corresponding counterterm 
diagrams. We depict the QCD box diagrams in Figj2] as a representative. The amplitude 
for the diagrams with gluon in loops may contain both UV and soft IR singularities, while 
the amplitude for the diagrams with gluino but no gluon in loops contains only the UV 
singularities. 

In order to remove the UV divergences, we employ the modified minimal subtraction 
(MS) scheme [16] to renormalize the strong coupling constant and the on-shell (OS) scheme 
to renormalize the relevant colored fields and their masses, respectively. After doing the 
renormalization procedure, the UV singularities are eliminated. Here we list the definitions 
of the counterterms of the wave functions of gluino, top squarks, top-quark, the mixing 
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Figure 2: The QCD box diagrams for the decay process g — > ttx\ ■ 

matrix elements of top squark sector, the strong coupling constant, the complex masses of 
top squarks and the mass of top quark as 
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where Uq and U* are the bare and renormalized 2x2 mixing matrices of the top squark 
sector, tifi and tj are the bare and renormalized fields of top squarks, and 5Z\- = 1,2) 
are the renormalization constants of top quark fields. 

In above counterterm definitions, we split the real bare mass of top squark squared (m| ) 
into complex renormalized mass and mass counterterm (i.e., fif and 5/i~). The bare fields 



of top squarks are split in complex field renormalization constants and renormalized fields. 
Similarly, U^ is separated into renormalization constant part and renormalized matrix 
element which can be complex. As a consequence, the renormalized Lagrangian, i.e. the 
Lagrangian in terms of renormalized fields and parameters without counterterms, is not 
hermitian, but the total Lagrangian keeps hermitian. 

B. Complex renormalization for top squark sector 

Because we use the CMS to deal with the possible top squark resonance, the normal 
one-loop integrals must be continued onto the complex plane. The formulas for calculating 
the IR-divergent integrals with complex internal masses in the dimensional regularization 
scheme are obtained by analytically continuing the expressions in Ref. [17] onto the complex 
plane. The numerical evaluations of IR-safe iV-point (N = 1, 2, 3, 4) integrals with complex 
masses, are implemented by using the expressions analytically continued onto the complex 
plane from those presented in Ref. [18]. In this way, we created our in- house subroutines 
to isolate analytically the IR singularities in integrals and calculate numerically one-loop 
integrals with complex masses based on the LoopTools-2.4 package [151 [19]. The analytic 
expression for scalar one-loop 4-point integral in the complex mass scheme is also given in 
Ref. [20J, we make a careful comparison between ours and that provided in Ref. [20]. It shows 
that they are in good agreement. 

In our calculation the external gluino and top quark are real particles, and their QCD 
on-shell self energies do not involve any absorptive parts. But they become complex via the 
complex top squark masses. The renormalization of top squark sector is more complicated. 
It involves the renormalization of the mixing of the two top squarks. 

With the renormalization conditions of the complex OS scheme [2~Tl |2"2"]. the renormal- 
ization constants of the complex masses and wave functions of top squarks are expressed 
as 



5/i|=4(/i|), 6Zl = -ZU4), 8Zl = J 3 _ * , (z,j = l,2,z^j). (3.7) 



We perform Taylor series expansions about real arguments for the top squark self-energy 





as 
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= 4(m|)-zm f T^s|(m|) + O (afj , (z,j = l,2), (3.8) 

where Sg(m?) = g, 2 | p2=m 2 , S^- ~ C(a s ) and (//|. - m|) ~ 0(a s ). By neglecting the 
higher order terms and using /i~ = m~ — im^. , we get approximately the mass and wave 
function renormalization counterterms of scalar top squarks as 

54 = E|(m|) + (4-m|)S|(m|), (3.9) 

Ml = -4K), (3.10) 
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By adopting the unitary condition for the bare and renormalized mixing matrices of top 
squark sector, Uq and U l , we obtain the expression for the counterterm of top squark mixing 
matrix 5U l as 



5U l = -(5Z i -8Z^)U*. (3.12) 
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With the definitions shown in Eqs.( l3.3| ) and f )3.5p . the counterterms of the top squark mixing 
matrix elements can be written in the following form: 

6U l = 7 E (*4 " E>iy (hi = 1,2). (3.13) 

4 k=i 

The explicit expressions for unrenormalizad self-energies of top squarks, E^(p 2 ) = 
1,2), are written as 
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From Eqs. fl3TT0|) . flSTTTj) . ( ED and fl3TT5|) . we can obtain Im(5Z^) = Im(<5%). The 
Eg. ( 13. 13]) tells us that we can choose both the renormalized mixing matrix elements {Ufj) 
and their counterterms {SU^) made up of real matrix elements. Then the matrix TP can be 
expressed explicitly as 

TP = ( C ° s6 t Sine A. (3.16) 
\ — sin u t cos u t I 



C. Renormalizations for g and a s 

Since in our LO and NLO calculations we adopt the complex mass scheme, the con- 
ventional top squark masses are replaced by the renormalized top squark complex masses 
fit. (i = 1, 2) everywhere, including in the expressions for the gluino self-energy and coun- 
terterm of the strong coupling constant. The renormalized one-particle irreducible two-point 
Green function of gluino is defined as follow 

if*(p) = i{£> - m s ) + i[frP L tl{ P 2 )+ i>P R ±* R tf) + P L £l L (p 2 ) + P R tl R (p% (3.17) 

where Pl,r = |(1 T 7s)) S|> ^r, £§l and E§ fl are renormalized gluino self-energies. By 
taking u-, d-, c-, s-, 6-quarks being massless, we express the corresponding unrenormalized 
gluino self-energies explicitly as, 
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By adopting the OS scheme to the external gluino, the wave function renormalization con- 
stants of gluino can be fixed by using the following formulas 
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For the renormalization of the strong coupling constant g s , we adopt the MS scheme 
at the renormalizatiion scale /z r , except that the divergences associated with the top quark 
and colored SUSY particle loops are subtracted at zero momentum [21]. We define the 
counterterm of the strong coupling constant 5g s as a summation of the SM-like QCD term and 
SUSY-QCD term (i.e, 5g s = §g( SM ~ hke ) _|_ 5g( SQCD ^), and these two terms can be expressed 
as 
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where the summation is taken over the squark indices of w, d, c, s and 6. The notations 
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The number of colors N = 3, the number of light flavors nif = 5, and 4 = — j E + ln(47r). 
There are two regularization schemes, dimensional regularization scheme and dimensional 
reduction regularization scheme, customarily used in a supersymmetric gauge theory. It 



is well-known the dimensional reduction regularization scheme preserves supersymmetry at 
least to one-loop order, therefore, this scheme is a proper regularization scheme for per- 
forming the renormalization in supersymmetry In this paper we adopt the dimensional 
regularization scheme to calculate the NLO QCD corrections. This scheme is easier than 
the dimensional reduction scheme to handle in general, but violates supersymmetry because 
the number of the gauge-boson and gaugino degrees of freedom are not equal in f ^ 4 di- 
mensions. To subtract the contributions of the false, non-supersymmetric degrees of freedom 
and restore supersymmetry at one-loop order, a shift between that the qqg Yukawa coupling 
g s and the qqg gauge coupling g s must be introduced [25] 




where Ca = N = 3 and Cp = 4/3. Similarly the electroweak t{tx\ couplings should take 
a finite shift also, which may be written as ae + bY t , with e = e[l — a s Ci?/87r] and Y t = 
Y t [l — 3a s CV/87r] in terms of the electromagnetic coupling e and the quark-Higgs Yukawa 
coupling Y t oc em t [25J. In our numerical calculations, we take all these shift into account 
and keep Tnlo only up to 0(a 2 s a ew ). 

For the renormalization of the wave function and mass of top quark by using OS scheme, 
we use the expressions Eqs.(10)-(13) in Ref.[26] with the replacement of m t ~. — > The 
0(a s ) QCD virtual corrections to the decay g — > ttxi in the MSSM can be expressed as 



where Aitree is the Born amplitude for the decay g — > ttx\ , and A4 ne-ioo P is the renormalized 
amplitude of all the one- loop level Feynman diagrams involving the virtual gluon or/and 
gluino. 

III. 2 Soft gluon emission corrections 

In our calculations we denote the real gluon emission decay process g — > tix\g as 
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Figure 3: The tree-level Feynman diagrams for the real gluon emission decay process g — > 

and the corresponding tree-level diagrams are shown in Figj3] 

The real gluon emission decay g — > tixlg contains only soft IR singularity, but no collinear 
IR singularity due to gluino, top quark and top squarks being massive. This singularity can 
be conveniently isolated by slicing the phase space into two different regions defined by an 
arbitrary small cutoff. This method of dealing with soft IR singularity is called phase space 
slicing (PSS) method [27]. By introducing a small cutoff 8 S , the phase space of g — > ttxlg is 
separated into two regions, according to whether the emitted gluon is soft, i.e., E 5 < 5 s m g /2, 
or hard, i.e., E 5 > 5 s rrig/2. Then the decay width of the real gluon emission process can be 
expressed as the summation of the contributions over the two phase space regions, i.e., 

AT real (g ->• ttxlg) = AT soft {g -> ttj&g) + AT hard (g -> ttftg), (3.28) 

where AT so f t is obtained by integrating over the soft region of the emitted gluon phase space 
and contains the soft IR singularity. AT hard is finite and can be evaluated by using Monte 
Carlo technique in four dimensions [28]. Although AT so f t and AT hard depend on the soft 
cutoff 8 S , the total decay width of the decay process g —> ttxlg , Ar rea /(^ — > ttx\ + g), 
is independent of the arbitrary cutoff S s . That independence is verified in our numerical 
calculations. In further calculations we set S s = 1 X I0" 4 . The decay width of the gluon 



bremsstrahlung g — > ttxig over the soft gluon region can be expressed as [29] 

Ar so/t = T LO <g> — 22 (Ti • T i >;./>/). (3.29) 

where Tj are the color operators [291 EHl El] , and gij are the soft integrals defined as 
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By using the definitions of color operators, we get the expression of T so f t as 
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Then the total decay width including the NLO QCD corrections of the decay process g — > 
ttxi can De obtained by summing all the contribution parts: 

"I - virtual + ^soft + ^hard- (3.32) 

IV Results and Discussions 

In this section we present and discuss the numerical results for the decay g — > itXi ■ we take 
two-loop running a s in both LO and NLO calculations. The number of active flavors Nf = 5, 
and the QCD parameter A^' IS = 226 MeV. We set the renormalization scale being equal 

TO-0 

to \i — no — mt H — ^ by default. We take the SM parameters as a ew {m 2 z ) = 1/127.916, 
m t = 172 GeV, m z = 91.1876 GeV and m w = 80.399 GeV [33]. 

Recently, the observations at the LHC indicate that the generic lower bound on gluino 
is 700 GeV at the 95% confidence level in simplified models containing only squarks of the 
first two generations, a gluino octet and a massless neutralino [31], and there are hints of 
the SM-like Higgs boson with rrih — 125 GeV [35]. If these hints for the Higgs boson mass 
are true, then that strongly suggests the top squarks are light. Accordingly, we consider the 
benchmark point SPSQ scenario, which is proposed in the SPA Convention and Project 
[TTJ [T21 [36J, as a numerical demonstration. The relevant masses of SUSY particles and 
parameters at the SPSQ point required in our numerical calculations are listed in Table [TJ 
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Table 1: Relevant SUSY parameters obtained by using IS A JET 7.82 [37] with the input 
parameters at the reference point SPS6. 

The parameters presented in Tableware the independent SUSY input parameters adopted 
in our calculations. All the other SUSY parameters needed in this work can be determined 
by those in Table [TJ For example, the elements of neutralino transformation matrix N can 
be given by M 2 , /i, tan /3 and m^o. Since we assume u-, d-, o, s-, 6-quark are all massless, the 
squark mixing matrices U q (q = u, d, c, s, b) should be all unit matrices. The elements of top 
squark mixing matrix V t are determined by % In order to investigate the dependence of the 
cross section on one of the parameters in Table [IJ we will vary the specific parameter while 
keep the other SUSY parameters in Table [1] fixed in the following calculations except for 
the top squark total decay widths. From our numerical calculations by using IS A JET 7.82 
with the input parameters at the SPS6 point [TTJ [12!, 36J , we get Y^im^ = 497.89 GeV) = 
2.98 GeV and Yi 2 {mi 2 = 678.64 GeV) = 9.70 GeV. For simplicity, we take the top squark 
decay width as T h {m h ) = m ^ GeV x 2.98 GeV {Y h {m h ) = 678 ^ 2 GeV x 9.70 GeV) when 
m h ( m t 2 ) ^ s variable. 

We show the LO, NLO QCD corrected g — > ttxl decay widths and the corresponding 
relative correction (5 = Vnl °~J lo ) as the functions of the renormalization scale in FigsJUa) 
and (b) separately, where we denote \i = fi r and /i = m t + m^o/2. If we define the scale 
uncertainty for the g — > tix\ process as r\ = l r (^- 5 w)~J^-°- 2 w)[ ^ f rom ^ e curves j n FigJU^a) 
we can figure out the corresponding uncertainties being r\ LO = 43.0% and r] NLO = 38.4%. 
It is clear that the LO decay width is strongly related to the renormalization scale in the 
plotted fi range, while the NLO QCD corrections obviously reduce the scale uncertainties. 
FigJUb) shows that the relative QCD correction S raises from —60.7% to —6.9% when the 



renormalization scale goes from 0.2/xo to 5/io, and the relative correction has the value of 
—28.2% at the location of fi = /x - 



0.20 




Figure 4: (a) The dependence of the LO and the QCD corrected g —> ttjli decay width on the 
renormalization scale, (b) The dependence of the relative correction for the g — > ttxl decay 
width on the renormalization scale. 

The decay widths of g — > tt\\ at 0(a s a ew ) and up to 0(a^a ew ), and the corresponding 
relative correction(£ = Tnl °~^ lo ) as functions of the ratio of the vacuum expectation values 
(VEVs) tan f3 are depicted in Figj5]^a) and Figl^b) separately, with tan (3 varying from 1 to 
35 and the other SUSY input parameters being adopted from the reference point SPS6 in 
Table [TJ The curves for the decay widths Tlo and Tnlo drawn in Figj5]^a) demonstrate that 
the NLO QCD correction suppresses the LO decay width of the decay process g —> ttxl ■ 
Both the decay widths Tlo and Tnlo quantitatively fall down when tan/3 goes up from 1 
to 5, and vary smoothly when tan > 5. Figj5]^b) shows the relative correction raises from 
—41.7% to —26.6% with the increment of tan/3 from 1 to 35. 

The dependences of the LO, NLO QCD corrected decay widths of g — > ttxi and the 
corresponding relative correction on the gluino mass m- g are depicted in Figsj6]^a-b), with 
the parameter m~ g in the range of [600, 1000] GeV and the other SUSY input parame- 
ters being adopted from the reference point SPS6 in Table [TJ Figj6](a) shows that the 
NLO QCD corrections are negative, and the LO and NLO QCD corrected decay widths 
at the position of m g = 721.8 GeV (1000 GeV) can reach 0.1490 GeV (1.364 GeV) and 
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Figure 5: (a) The LO and the QCD NLO corrected decay widths of g — > tixi as the 
functions of tan/3, (b) The corresponding relative correction as the function of tan j3. 




I ......... i ... I -0.4U ' ■ ' ' ' ■ ' ■ ' ' ' ' ' ' ' ' ' 

600 650 700 750 800 850 900 950 1000 600 650 700 750 800 850 900 950 1000 

m § (GeV) m § (GeV) 



Figure 6: (a) The LO and the QCD NLO corrected decay widths of g — > ttxi as the 
functions of gluino mass m g . (b) The corresponding relative correction as the function of 
m~ g . 

0.1069 GeV (1.103 GeV), respectively. In Figj6]Ja) there exist two knees on the NLO 
curve at the positions satisfying the conditions of m g = mf + m t = 669.9 GeV and 
m-g = mi 2 + m t = 850.6 GeV, but this resonance effect is not observable on the LO 
curve. Figj6]^b) shows that the relative corrections at the points of m g = 721.8 GeV 
and m g = 1000 GeV are —28.2% and —19.2%, respectively. We can see obviously from 
Figj6]^b) that the curve for the relative correction has two obvious spikes at the vicinities of 
rrig = 699.9 GeV and m g = 850.6 GeV due to the top squark resonance effects. 

In FigfTla) we plot the LO and the NLO QCD corrected decay widths, Ylo and Tnlo, as 
the functions of t\ mass by keeping the other SUSY input parameters at the reference point 




Figure 7: (a) The LO and the QCD NLO corrected decay widths as the functions of m^. 
(b) The corresponding relative correction as the function of m^. 

SPS6 shown in Table [1] except the varying with m% . The corresponding relative NLO 
QCD correction (5) as the function of is shown in FigJT^b). FigfT^a) shows that the curves 
for Ylo and Tnlo are significantly enhanced in the m% x range of [362.4, 549.8] GeV due to 
the resonance effect of the intermediate on-shell t\ in the decay chains of g — > tt\ — > itxi and 
g — > tt\ — > tix\i where the two turning points satisfy the conditions of = m^o + m t = 
362.4 GeV and = nig — m t = 549.8 GeV. The LO and the NLO QCD corrected decay 
widths reach their maxima of 0.2341 GeV and 0.1760 GeV at the position of mi = 440 GeV, 
respectively. On the curve of the relative correction in FigJT^b) there exist two spikes at the 
vicinities of m t - = 362.4 GeV and m^ = 549.8 GeV, which also reflects the resonance of the 
intermediate on-shell t% in the decay process g — > tix\- 

The LO, NLO QCD corrected decay widths and the corresponding relative NLO QCD 
correction as the functions of m^ 2 are depicted in FigJHl^a) and FigJH^b), respectively, with 
mi 2 varying in the range of [520, 700] GeV . In Figs|S^a,b) all the SUSY input parameters 
are taken from the reference point SPS6 except for m 4 ~ 2 and T^. FigJSJ^a) shows that the 
LO and the NLO QCD corrected decay widths have the values of 0.1586 GeV (0.1490 GeV) 
and 0.1204 GeV (0.1070 GeV) at m h = 520 GeV (700 GeV), respectively. We can see 
from this figure that in the region of m 4 ~ 2 < 549.8 GeV both Y LO and T NLO have relative 
large values due to the resonance effect of the intermediate on-shell t 2 (where the condition 



m t 2 — m g ~ m t = 549.8 GeV is satisfied). From FigJH^b) we can read out that the relative 
corrections are —24.1% and —28.2% at the positions of m 4 ~ 2 = 520 GeV and m ( - 2 = 700 GeV, 
respectively. Again, there is a spike on the curve for relative correction at the vicinity of 
mi 2 = 549.8 GeV induced by the resonance effect of the intermediate on-shell £2- 
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Figure 8: (a) The LO and the NLO QCD corrected decay widths as the functions of m t ~ 2 . 
(b) The corresponding relative correction as the function of m t - 2 . 

The decay widths of g — > tix\ at the LO and the QCD NLO, and the corresponding 
relative correction versus m^o are shown in FigfJ^a) and Figj9]^b) separately, with m^o being 
in the range from 100 GeV to 200 GeV and the other SUSY parameters being taken from the 
reference point SPS6 shown in Tabled] The curves for the LO and the NLO QCD corrected 
decay widths decrease with the increment of m^o. The curve for Ylo (Xnlo) goes down 
from 0.211 GeV (0.152 GeV) to 0.139 GeV (0.0964 GeV) when m^o varies from 100 GeV 
to 200 GeV. The corresponding relative corrections at the points of m^o = 100 GeV and 
m^o = 200 GeV are -28.2% and -30.9%, respectively. 

The LO and the QCD NLO corrected distributions of the top-pair (tt) invariant mass M t i 
and the missing energy E miss (^, ^ff, -§i& and g^ft) at the reference point SPS6 
are shown in Fig JTOT a) and Fig JTOT b). respectively. The kinematics of the decay process 
g —?• tixi at this point constraints the dr ^ JV£CI distributions in the range between M™ 71 = 



2m, ~ 340 GeV and ML 



530 GeV as shown in Fig JTOT a). and the 
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distributions are limited in the range between E 7 ™^ = m^o ~ 190 GeV and E, 
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Figure 9: (a) The LO and the QCD NLO corrected decay widths of g — > ttxl as the 
functions of m^o. (b) The corresponding relative correction as the function of m^o. 

(4m2- m 2 ) 

— 2m _ 1 ~ 304 GeV as shown in Fig lTOT b). We compare these differential decay widths 
and find that the NLO QCD corrections suppress the corresponding LO differential decay 
widths and significantly, but do not change obviously the line shapes of the LO 

differential decay widths. From Fig JTOT a) we see that the M t i distributes mainly in the range 
of [350, 480] GeV, and the differential decay widths raise slowly with the increment of M t i 
in this region. Fig lTOT b) shows that the missing energy events are mainly concentrated and 
nearly uniformly distributed in the E miss range of [225, 300] GeV for both the LO and QCD 
NLO corrected distributions. 



V Summary 

In this paper we calculate the NLO QCD corrections to the decay process g — > tixi m t ne 
MSSM. As a numerical demonstration we present and discuss the NLO QCD corrections 
around the SPS6 benchmark point. We find that the NLO QCD corrections significantly 
suppress the corresponding LO decay width, and at the scenario SPS6 point the LO and the 
NLO QCD corrected decay widths have the values of 0.1490 GeV and 0.1069 GeV respec- 
tively, and the corresponding relative correction is —28.2%. We analyze the dependence of 
the NLO QCD corrected decay width and the corresponding relative correction on the ratio 
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Figure 10: (a) The tt invariant mass distributions at the LO and the QCD NLO ( dj ^° , 
dT dM L ° ) ^ or ^ ne SPS6 parameter set in Table [TJ (b) The missing energy distributions at the 
LO and the QCD NLO {j§Mk, with the SPS6 parameter set listed in Tabled! 

of the vacuum expectation values, gluino mass, scalar top-quark masses and the lightest 
neutralino mass, respectively, around the scenario point SPS6. Our numerical results show 
that the absolute NLO QCD relative correction can exceed 30% in our chosen parameter 
space. Therefore, it is necessary to take the NLO QCD corrections into account for the 
precise experimental measurement at future colliders. We compare the distributions of the 
tt invariant mass and the missing energy at the LO and the QCD NLO and find that the 
line shapes of the differential decay widths at the LO, and f^? ss , are not obviously 
changed by the NLO QCD corrections. 
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